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Abstract. We present a first step toward a microscopic theory for the glass transition 
in systems with trivial static correlations. As an example we have chosen N infinitely 
thin hard rods with length L, fixed with their centers on a periodic lattice with lattice 
constant a. Starting from the iV-rod Smoluchowski equation we derive a coupled 
set of equations for fluctuations of reduced fc-rod densities. We approximate the 
influence of the surrounding rods onto the dynamics of a pair of rods by introduction 
of an effective rotational diffusion tensor D(r2i,r22) and in this way we obtain a self- 
consistent equation for D. This equation exhibits a feedback mechanism leading to a 
slowing down of the relaxation. It involves as an input the Laplace transform vo{l/r) at 
z = 0, Z = i/a, of a torque-torque correlator of an isolated pair of rods with distance 
R = ar. Our equation predicts the existence of a continuous ergodicity-breaking 
transition at a critical length Ic = Lc/a. To estimate the critical length we perform 
an approximate analytical calculation of vo{l/r) based on a variational approach and 
obtain Z™'' ~ 5.68, 4.84 and 3.96 for an sc, bcc and fee lattice. We also evaluate 
vo{l/r) numerically exactly from a two-rod simulation. The latter calculation leads to 
^num ^ 3.45^ 2.78 and 2.20 for the corresponding lattices. Close to Ic the rotational 
diffusion constant decreases as D{1) ~ {Ic — l)'^ with 7=1 and a diverging time scale 
tt ~ \lc — , (5 = 2, appears. On this time scale the t- and /-dependence of the 
1-rod density is determined by a master function depending only owt/t^. In contrast 
to present microscopic theories our approach predicts a glass transition despite the 
absence of any static correlations. 
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1. INTRODUCTION 

In the last two decades huge theoretical efforts were made in order to describe the glass 
transition on a microscopic level in systems without quenched disorder. Microscopic 
means that, e.g., the glass transition temperature and the dynamical properties in 
its vicinity can be calculated from equations which use as an input the knowledge of 
the interactions between microscopic constituents of the system (atoms or molecules; 
hereafter we use the term particles). Most of this activity has been devoted to structural 
glasses |jl], ||, ^, ||. Although experiments f^, 0, || have demonstrated that plastic 
crystals can exhibit an orientational glass transition which has many features in common 
with the structural glass transition there seems to be no analytical theory for such 
systems. |f] On the other hand Molecular Dynamics (MD) simulations [0, [12|, |1^ 



have confirmed the similarity between glassy behavior of plastic crystals and that of 
supercooled liquids. 

The first successful microscopic theoretical approach to glassy dynamics in 
supercooled liquids came from so-called mode coupling theory (MCT) ^ |[ . This theory 
was first suggested and mainly worked out in details by Gotze and his coworkers. For 



reviews the reader may consult refs. |]T3|, 0, |T6[ . MCT derives a closed set of equations 



for time-dependent correlators, like the intermediate scattering function S{q, t). This set 
uses as an input the corresponding static correlations, e.g. S{q), which depend only on 
thermodynamic variables like temperature T, density n, etc. The static correlations can 
be calculated from the microscopic interactions by use of either analytical or numerical 
tools |T^. By decreasing T (or increasing n) the peaks of S{q) grow. If their magnitude 
becomes large enough a transition from an ergodic to a non-ergodic phase occurs at a 
critical temperature Tc (or density Uc). This dynamic transition is interpreted as an ideal 
glass transition. MCT makes several nontrivial predictions for, e.g., time dependence 



of correlators |jTj, |T5|, [T^ which have been successfully tested (see |]T8|, |T9], [2^, ^ for 
reviews) . 

The second microscopic theory of the glass transition has been developed by Mezard 
and Parisi 0. Their first principle theory combines the use of the replica idea known 
from the spin glass theory and of the liquid state theory techniques. It is related to 
an earlier density functional approach [Q. Mezard and Parisi's theory starts from the 
microscopic Hamiltonian for m identical replicas of the system to which a symmetry 
breaking field is added. The free energy of the replicated system is then calculated 
and the configurational entropy of the original system is extracted from it. The theory 
predicts the existence of a critical temperature Tf (which is below Tc of MCT) at which 
the configurational entropy vanishes. This is then interpreted as an ideal static glass 
transition. 

Besides these two microscopic theories, let us mention one other attempt to describe 
the glass transition that does not resort to microscopic potentials: based on several 

f The microscopic theory developed by Michel and coworkers [T^ is valid for the orientational glass 
transition of mixed crystals, i.e. crystals with quenched disorder. 
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qualitative assumptions Edwards and Vilgis |]22| obtained the dependence of the self- 
diffusion constant D for a system of hard rods on physical control parameter x, like, 
e.g., concentration, diameter, or length of the rods. A dynamical glass transition occurs 
when D becomes zero. A mean field treatment results in D ~ Xc — x and an improved 
approach leads to the Vogel-Fulcher law D ~ exp[74/(xc — x)]. Close to the critical 
point Xc a power law dependence of D on [xc — x) is found with a non-trivial exponent 
7 = 7/6. 

Both above described microscopic theories have at least one feature in common: 
non-trivial equilibrium correlations. These correlations are of primary importance for 
both theories. For example, MCT requires that the vertices which enter into the self- 
consistent equations for the correlation functions must depend on at least one control 
parameter. Recently, it has been shown that a colloidal suspension with hard-core plus 
attractive Yukawa potential can undergo a gel transition in a special limit: — > 
and K ^ oo such that F = K'^ip/h = const. Here is the volume fraction, K the 
interaction strength and b a measure of the inverse range of interaction. In that limit 
the direct correlation function c{q) is proportional to K and becomes independent of 



[^. This implies that the static correlations vanish, i.e. S{q) 1. But since 
the vertices are proportional to and involve a bilinear product of c{q) they become 
proportional to F. Accordingly, despite the vanishing of the static correlations there is 
still an ergodicity breaking transition driven by the increase of F |2^, ^ Of course, 
this type of behavior is not generic in the (p — K phase space. In addition, whether the 
mode coupling approximation remains reasonable must be still checked. 

One may now ask the following question: Can there be a generic glass transition if 
and only if there exist nontrivial static correlations? To find an answer let us consider 



a liquid of hard rods. This system was used as a paradigm in ref. E2]. If the rods 
thickness d is finite the static correlations are non-trivial and a glass transition will be 
predicted by either theory. However, in the limit of infinitely thin hard rods, i. e. d = 0, 
with finite length L and finite concentration c, the static properties become trivial. Note 
that in this limit there is no transition to a nematic phase. Also, according to either 
of the two microscopic theories there is no glass transition! In the case of MCT the 
vertices vanish. The absence of a glass transition is also consistent with a microscopic 
theory worked out by one of the authors and his coworkers |26, 27 . There, it has been 
shown that for infinitely thin hard rods with randomly frozen orientations translational 
diffusion never ceases if the diffusion constant along the rod's direction is nonzero. 

Now suppose we fix the infinitely thin hard rods with their centers on the sites of 
a three-dimensional periodic lattice. Such a model has already been suggested earlier 



2^ , |28[] . It can be used to describe the steric hindrance of the orientational degrees of 
freedom in plastic crystals [^, 0, ^ |TT|, . Again, both MCT and the replica theory 

exclude a glass transition, due to the lack of static correlations. The fact that MCT 



in its present form does not predict a glass transition was already stressed in ref. 



Nevertheless, computer simulations PO, RDl showed glassy behavior and suggested a 
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dynamical ergodi city-breaking transition at a critical length L^. || This simple example 
makes clear that non-trivial static behavior is not necessary for a glass transition. It is 
the main goal of our paper to present a microscopic approach to the glass transition for 
systems without static correlations. 

A short account of our approach for the case of a simple cubic lattice was given in 
ref. |^T|. In the present paper we will discuss our model and the theoretical framework 
in more detail and will extend the results to bcc and fee lattices. Furthermore, we will 
investigate the type of glass transition and the long time dynamics in its vicinity. 

The outline of our paper is as follows. In the next section we will present and 
discuss the model. The third section contains the analytical theory used to describe the 
glass transition. The results are given in the fourth section and the final section includes 
a summary and some conclusions. 

2. MODEL 

We consider N hard, infinitely thin rods [i.e. diameter d = 0) with length L and with 
their centers fixed at the lattice sites of a periodic lattice with lattice constant a. The 
dimensionless control parameter is the reduced rod length, I = L/a. In the following we 
will restrict ourselves to a simple cubic lattice. Its unit cell is depicted in Fig. |I|. All 
qualitative conclusions which will be drawn in this paper remain valid for all periodic 
lattices. 

It is obvious that for / < 1 each rod can rotate freely: there is no steric hindrance. 
For / above one, collisions between a rod and adjacent ones appear. If 1 < / < -\/3 
a given rod interacts only with the six nearest neighbors. Further increase of / will 
introduce collisions between next-, third-, fourth-, etc. nearest neighbors, generating 
a strong increase of steric hindrance. Since we have chosen infinitely thin rods there 
is no equilibrium transition to an orientationally ordered phase. On the other hand, 
due to the rise of steric hindrance with / we expect an orientational glass transition as 
detected in ref. and Although steric hindrance already exists for /> 1, it 

seems that no ergodicity-breaking transition can occur for / < 2. Let us consider the 
following simple argument. If Z > 2 then the nearest neighbor rods of a tagged rod 
(rod 1) can "overlap" with, e.g., the top of the tagged rod (see Fig. 2). This induces a 
"cage" for the tagged rod. Note that this argument is not rigorous, because for / < 2 
there might exist other configurations similar to the one shown in Fig. ^ which could 
block the tagged rod's rotation. But the value / = 2 also plays a role from a different 
point of view. Consider two neighboring rods with their centers along the z-axis. Their 
orientations can be described by Vti = {6i,(j)i) and Qj = {6j,(f)j). At a collision it is 
either (f)j = (f)f (see Fig. ^ a) or (pj = 0f + vr (see Fig. ^ b), for fixed 6*^, 9j and 0,. The 
notation (pf and 07 means that rod j is behind and in front of rod i, respectively (see 
Figure Now, it is easy to prove that for arbitrary 6i and 6j there exists only one 

I Not that different cubic fattices were used in refs. |^ and |3^; moreover tfie autiiors of ref. | ]30| 
fixed tlie rods witli one of their endpoints at the lattice sites. 
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contact (either at (j)f or at (j)f + vr) for / < 2 (see Fig. |^ a,b). In that case can be 
varied freely by 27r. But for I > 2 there is a finite range for both angles 6i and 6j such 
that two contacts (at (pf and (j)f + vr) exist (see Fig. |^ c). Therefore the "phase space" 
[0, 27r] for — is decomposed into two "ergodic" components [0, vr] and [vr, 27r] for 9i 
and 6j kept fixed. Transitions between both components are forbidden. If the rods were 
fixed at the lattice sites with one of their end points (like in ref. then there is only 
one contact, for all I > 1. 

The condition for a collision can be quantified as follows ||3^. Let Uj and Uj be the 



unit vectors along rod i and j and r^j the vector connecting the center of rod j with 
that of rod i (see Fig. Uj and Uj define a plane. Let be the component of Tij 
perpendicular to that plane and rf^ = |r^|. Denote the distance of the point of contact 
to the center of rod i and j by, respectively, Sij and Sji (Fig. §). Then a collision occurs 
if: 



r ■ ■ 

1-] 



0+ , \si.j\ < L/2 , \sj,\ < L/2 (1) 



where 0"*" stands for limit d — 0. Simple algebra yields: 

Sij = -[(uj • Tij) - (uj ■ Uj){uj ■ rij)]/[l - (uj ■ Ujf] (2) 

Sji = [{Uj ■ Vij) - (uj ■ Uj){ui ■ rij)]/[l - (ui ■ UjY] (3) 

Note (i) that Sij and sji can be negative; in that case the contact points into the direction 
of — Uj and —Uj, respectively and (ii) that Vij is time-independent and is always equal 
to a lattice vector R^- = Rj — Uj. 

The configurational part of the partition function is given by 

N 

Z,{T, N;l)= Y[ dn, e-/^^^^--'^-) = {Anf , /3 = (fc^T)"! (4) 
i=i 

since the hard rod interaction potential V{Qi, ■ ■ ■ , ^at) vanishes almost everywhere. The 
corresponding contribution —/c^TiV In 47r to the free energy is analytic in temperature 
T and does not depend on /. Hence, there is no equilibrium phase transition. The 
equilibrium pair distribution function gij\^li, Q.j) is related to the probability to find a 
rod at Rj and Rj(7^ Rj) with orientations Qi and flj, respectively. It is given by: 

^If (r^., Q,) = 1 - ^ (o+ - ri^) e (I - e - (5) 

where the Heaviside functions are nonzero at a collision, which happens if conditions 

(9) 

Eq. ([1|) are fulfilled. Since a collision is non-generic, glj (^li,nj) = 1 for almost all fij 
and Qj. Hence there are no static correlations, except on a set of measure zero. 
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In this section we will describe the dynamics for our model and the approximations 
leading to a closed set of equations. We assume that the microscopic A^-rod dynamics 
of the system is Brownian rather than Newtonian. Previous studies of systems with 
nontrivial static correlations indicated that glassy dynamics should not depend on that 



choice H, 

The starting point of the theory is the so-called generalized Smoluchowski equation 
for the A^-rod probability density P/v(^i; ' ' ' i^N\t) 



d ^ 

— P^(fii,---,l]^;t) = Do E V„-[V„-T„(f]i,---,l]^)]Pjv(^^i,---,f^Jv;t)(6) 

^'^ n=l 

where Dq is the bare rotational diffusion constant, V„ = Vn„ = u„ x Vu„, and 

T„(fil, ■ ■ ■ , ^n) = E Tjn(^i) ^n)- (7) 



T^j(^]^,^7j) describes the singular torque when two rods at sites i and j with vector 
distance Yij collide: 

T,,-(n,,fi,) = .,,(u, X vj^)5{ri^ - 0+)^ (I - \s,,\) 6 - |.,,|) (8) 

where f^- = rjj/r^. Using Eq. (||) it is easy to prove that 

Vi^if {n.,Qj) = T,,{n„Q,)gif{Q„n,). (9) 

Eq. (H) describes the noncrossability condition. Note that Tij{Qi,Qj) differs slightly 
from Tij{fli,flj) used in ref. since there V„. was used instead of Vq.. Also, a 
similar generalized Smoluchowski equation had been used before to describe a liquid of 



infinitely thin hard rods with randomly frozen orientations p6|, ^ 



In a next step we present equations of motion for the reduced fc-rod density 



rn 



= J H rf^]„p^(fil,■■■,^]^;^). (10) 



These equations form an infinite, coupled hierarchy ||3^ of the following form 



-^0 E E v„ ■ 1 dQk+i Tn7ife+i (^rt) ^fc+l) 

xpl\~^Xjn„.--,n,^i;t) (11) 

with the k-Tod Smoluchowski operator, 

£('=) = DoE V„-[V„- E T„„,(fi„,l]„0]- (12) 
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and Ik = {rii, ■ ■ ■ , n^}. As initial conditions we choose: 

pS..„^ (l^i, ■ ■ ■ , fi,; 0) = 6{n,\no) gitnS^u ■■■,^k), (13) 

where S{fl\Q') = sm96{9 - 6') 6{(l) - 0'). Here, g^_^...^^{Qi, ■ ■ ■ , Qk) is the equihbrium 
/c-rod distribution function which is equal to one almost everywhere. 

To proceed we introduce the fluctuations of the k-Tod density for k > 2 



5p^S-n,{^i^ ■■■,^k;t)= pS..„,(fii, ■■■,Qk;t) 

-T-^TT gi'l-n^i^u ■■■,^k) pS(^^i; t). (14) 
(4vr) 

These fluctuations vanish in equilibrium and for t = 0. 

Substituting p^^^ from Eq. ([T^) into Eq. (|llD and taking into account that the 
equilibrium k-Tod rotational current density vanishes we obtain the following equation 

|pff (fii;t) + Vi-j«(fii;t) = o (15) 

where the 1-rod rotational current density reads 

}S{^i;t) = -Do{v,pS{n,;t)- 

Y: J dn2Tn,n,in,,Q2)SpSn,i^l,^2;t)\ . (16) 

At the next level, A; = 2, it follows by use of Eq. 

+Do{Vi ■ [Vi - T„,„,(fii,fi2)] + (1 2)}5pg)„^(^]l,^]2;t) + 

+A(^U^„n,;t) (17) 

where 



A^S^{Q^,Q2;t) 



gfl^{n,,n,)Y: ^ / rfn3T„,„3(f^i,^3)^pSl,(^^i,^^3;t) 



:i8) 



Let us interpret these equations: Eq. (|I5|) is the continuity equation. The 
corresponding 1-rod current density (Eq. (p!6|)) consists of two parts. The first describes 
the contribution from the free Brownian dynamics of a tagged rod at site rii and the 



GLASS TRANSITION IN SYSTEMS WITHOUT STATIC CORRELATIONS 



8 



second is due to the interaction of the tagged rod with a second rod at side n2- The 
first hne on the r.h.s. of Eq. (0) originates from the last term on the r.h.s. of Eq. (p^), 
its second hne describes Brownian dynamics of an isolated pair of rods at sites rii and 
n2, and the last line contains the influence of a third rod at 713 on the rods at ni and 
n2. In order to close this set of equations we follow the strategy of refs. and 
approximate the influence of a third rod by introducing an effective, nonlocal (in angular 
space and time) diffusion tensor T)^^{yt^VL']t) and simultaneously neglecting A^^^^ 



in2 • 



{DoVi ■ [Vi - T„,„,(fii,fi2)] + (1 2)}5pg)^^(^]l,fi2;^) ^ 

Vi ■ Jdt' I dn[B''''{n,, n[; t-t')- [v; - T„,„,(fi;, fi2)]5pSU(^^'i, ^^2; t)+ 





t 



+V2 ■ j dt' j rffip'=f^(i72, t-t')- [v'2 - T„,„,(fi^, n,)]5p^^l^{n^, a^- 1) 


= : (A(2)*5p(2))„^„^(f],,(]2;t), (19) 

A(%^(fi„fi2;t)^0. (20) 

The generalized rotational diffusion tensor D(f2, t) is defined through a 
constitutive equation relating the 1-rod current density and the (angular) gradient of 
the 1-rod density 

t 

j(\)(fi; t) = -Jdt'J dQ'B{Q, Q'; t-t')- Vn'pSi^'; t)- (21) 



As the final approximation we impose a self-consistency condition, 

D"^ = D. (22) 

Eqs. (p^SHT^), (pJH2^) form a closed set of equations for p^^\ Sp^"^^ and D. Taking 
their Laplace transform ^ it is straightforward to (formally) eliminate p'-^^ and 5p^^^ 
which results in: 

Do^kii^hz) 

+ E^/ c^^^3T.,„3(^i,^3)[(;^-A^'))-'*(T.,„3^7il-jS)kn3(f^i,^ 

= I dn, D-i(fii,fi3;^) -InA^s-.z) (23) 



where ^ denotes the Laplace transformed quantities. From Eq. (^) we find immediately 
the self-consistent equation for D(fii, Q2', z): 

§ Here we use the definition f{z) = J dt f [t) zt) , Re z > and take 5pnin.2{^i,^2\t = 0) = 



into account. 
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+ i E Tn,n3 (^^1, ^^2) [{Z - * (n,„3 g'i^l,)Un, {^U ^2, z) (24) 

Taking the limit z ^ in Eq. (^4]) and operating with (l/47r) / dVli J dQ2 on both 
sides of the resulting equation yields the following for the diffusion tensor D(f2i, ^2) = 
B{Qi,Q2;z = 0): 

47r(D-i) = Do 1 - * T)„,„3 TlJ (25) 



1713/ 

"3 



with 



{fnin2 hn-in2) 

= (4^ Z"^^^/ ^^2^7£U(^^l,1^2)/:,„,(l^l,^^2)Vn.(^^l,1^2), (26) 



and 



(Do)°^ =Do5"^. (27) 
T* is the transpose of T and for the adjoint 2-rod operator we find from Eq. (p!9|): 



(A(2)t*/)„^„^(fii,^]2) 

= [Vl + TnmA^lM] ■ J dns^{ns,^l) ■ V3/„,n,(f^3,^^2) 

+ [V2 + Tn,nA^2,ni)]- 1 rffigD (1^3, ^^2) " Vs/n.n, (^1 , ^s) • (28) 

The reader should note that (i) 5f^^^^(f2i, ^2) in Eq. (^61) can be skipped, because it is 
equal to one almost everywhere; of course, it must not be dropped if Vi or V2 act on it, 
like in Eq. ([T8|); in that case one can use Eq. (||); (ii) (D~^) does not depend on rii and 
n2 and (iii) the tensor (((A'^^^"'')"^ * T)„j„3 T^^^^) in Eq. (p5|) depends only on R^g — 
because of the lattice translational invariance. Accordingly " ") does not depend on 

rii. Eq. (p5D can also be rewritten as follows 

00 

47r(D-i) = Do ^ + E / dt{{e^''''' * T).,„3 T^^J. (29) 







Let us discuss Eq. (ppj). It is a functional equation for the rotational diffusion tensor 
D, since A'-^-*^ also involves D. It / < 1 then T„j„2(i7i, ^2) = and therefore Eq. (^) 
implies that 

D(^]l,^]2) = Do5(fii|i^2) (30) 

as it should be. Increasing / beyond one leads to an "increase" of the friction 
tensor (D~^), due to the positive definite, time-dependent torque-torque correlation 
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tensor ((exp(A(^)"''t) * T)„^„^ '^ning) which determines the second term in Eq. (|29|). 
This "increase" of the "renormahzed" friction tensor imphes a decrease of the 
"renormahzed" diffusion tensor D. Since this one enters into the exponent of exp 
the relaxation of the torque-torque correlation will slow down. This in turn leads to 
an "increase" of the second term of Eq. ( p9D leading to a further "increase" of 
and so on. This feedback mechanism, which is different from but still resembles that of 



MCT [|1J, [1^, |16|, finally may lead to the vanishing of diffusion and therefore to a glass 
transition. 

Although the functional equation possesses a rather clear structure it probably 
cannot be solved exactly. For this one would have to determine the eigenvalues and 
eigenfunctions of A^^^'''. Due to the singular character of T„j„2(fii,fi2) (cf- Eq. (§)) on 
the three-dimensional contact-hyper surf ace defined by Eq. (|l]) this does not seem to be 
feasible. In order to make progress and to get explicit results for D we will use the 
following additional approximation: 

(0(^1,^2))"^ ~ D(0rt(^^i|^^2). (31) 
Substituting Eq. (pTD into Eq. (^) and acting with | I] on both sides yields: 

D{l) = D,[l-v{l)] . (32) 

The 1-dependent coupling function v{l) follows from the second term of Eq. ( [29| ) by use 
of Eqs. (HD and (H): 



1 °° 

^(0 = ttE /rft(T.,.3-e^''''*T„.n3) (33) 



with the adjoint of the 2-rod Smoluchowski operator (cf. Eq. (|T2|)): 

£(2)t = + T„,„,(fil,^]2)] ■ Vi + (1 ^ 2) (34) 

where Dq is replaced by one. The calculation of the relevant quantity v{l) and the 
determination of a critical length 1^ at which D{1) vanishes will be presented in the next 
section. 



4. RESULTS 



Glass transition singularity 

Using the local approximation Eq. (|31|), we obtained Eq. (|3^ ) which constitutes a rather 
simple result. Before we calculate v{l) let us discuss D{1) and v{l) on a qualitative level. 
Since T„j„2(^i)^2) = for / < 1 it follows that v{l) = and therefore D{1) = Dq for 
/ < 1, as it should be. Next, we introduce: 

CO 

voil/rn) = \j rft(Ton • e^''^'* Ton) (35) 
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where r„ = |R„|/a. This is related to the torque-torque correlator of one rod at the 
origin and another one at site n. Note that the r.h.s. of Eq. ( pS]) depends on |R„|, but 
not on the direction of R^. vo{x), of course, vanishes for x < 1, is positive for x > 1 (see 
below) and we expect it to converge to a finite limit v"^ for x — oo. is completely 
determined by vq{x)\ 

v{l) = J2 Ml/rn) (36) 

n 

where the sum is restricted to such n for which Vn < l- This relation can be used to 
determine the asymptotic behavior of v{l) for I —>■ oo. Approximating the sum by an 
integral which becomes more and more accurate with increasing r„ we get: 

v{l) ^ J d^rvo{l/r) = An J dr r'^Vo{l/r) 

l<r<l 1 



An(^l^dxx^vo g)j 



I' 



which gives 



with 



v{l) ^c-l^ + 0{f) , / ^ oo (37) 



1 



47r / dxx^ t;o ( - ) > . (38) 







X 



That v{l) ~ for / — oo is obvious since all rods within a sphere of radius of 
order I will collide with the central one. Their number, of course, is proportional to l^. 
Because v{l) = for / < 1 and v{l) increases linearly with for large / there must exist a 
critical length Ic for which v{lc) = 1 and therefore D{lc) = 0. To determine Ic at which a 
dynamical glass transition takes place we have to calculate Vo{x). Although this quantity 
is much simpler than the second term in Eq. (^9]) , it cannot be determined exactly 
analytically (see below). Therefore the exact evaluation can only be done numerically. 
For such a numerically exact calculation we performed a simulation of the Brownian 
dynamics (defined by /Z^^-*^) of an isolated pair of rods with distance r„. The time- 
dependent torque-torque correlator resulting from this simulation is shown in Fig. |^ for 
several values of //r„. For t — > we find a power law divergence as expected and 
already proven for the force-force correlator of hard spheres [^, 0. Since only two 
rods are considered, the system is ergodic and therefore the torque-correlations relax 
to zero. Figure ^ shows that the relaxational behavior changes qualitatively around 
l/vn = 2 from a fast decay like for free Brownian dynamics to a rather slow relaxation. 
This crossover at //r„ ^ 2 is related to the properties discussed in section 2. Whether 
the long time decay of the torque-torque correlation is purely exponential, proportional 
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to t " exp(— At) or even a single power law t " (as it is true for two hard spheres p6| ) 



cannot be decided from our numerical data. Although the log-linear representation in 
Figure 5 exhibits a bending for large t (which would indicate deviation from a pure 
exponential) the exponential behavior could appear on a much larger time scale on 
which the statistical fluctuations of the numerical data prevent the determination of the 
precise long time decay. The numerical evaluation of the time integral in Eq. ( P3| ) leads 
to Vq^^{x) presented in Figure |^. The reader should note the strong increase starting 
at X = //r„ ^ 2. 

Now let us turn to the analytical calculation of Vq{x). For this we rewrite Vq{x) as 

i;o(x) = ^(Ton-(-/:(')^)-'Ton). (39) 

which is (up to the factor 1/3) the Laplace transformed torque-torque correlator at 
z = 0. 

Analogous to refs. ||26|, ^ we introduce a vector function f„^„2(^i) ^2) such that: 
£(2)tf^,„,(fii,l]2) = T„,„,(fii,fi2). (40) 



Eq. (pOD can be decomposed into a regular part 

[V? + V^]f„,„,(fii,fi2) = 0, (41) 
and a singular one: 

[Si2(ui X r^2) ■ Vl + S21(U2 X f^J ■ V2] f„i„2(r]l,fi2) = Si2(ui X V^^).{A:2) 



Eq. ( ^21) is to be satisfied for VLi,VL2 located on the three-dimensional contact 
hypersurface. Accordingly, Eqs. ( p] - P^ ) describe a boundary value problem with 
boundary value on the three-dimensional hypersurface embedded in the four- dimensional 
space built by the surface of two unit spheres. This is a difficult mathematical problem 
which seems to resist a rigorous solution. If it could be solved then the calculation of 
Vq{x) is reduced to the calculation of an integral over Vti and ^2: 

t;o(a;) = -i(Ton-fon). (43) 

Therefore, we resort to an alternative. It is easy to prove (by the use of Eq. (^)) that 
Eq. ( ^0|) is the variational equation of the functional: 



For a similar discussion see ref. It can be easily shown by partial integration that: 



/f . £(2)tf \ ^ 

+ {'^2f:,n,i^U^2)r] (45) 



GLASS TRANSITION IN SYSTEMS WITHOUT STATIC CORRELATIONS 



13 



which is negative for all f„j„2 7^ const, i.e. is a negative operator. Therefore 

one can follow ref. and prove that the solution i^^^^{yLi,VL2) of Eq- (iO) 



non-degenerate maximum of JF. From Eqs. (^), (^) and ( PD we then find that 

t;o(x)=^Kf„,(fii,r]2)], (46) 

i.e. the maximum value of JF is just 1^0(2^)- Since the maximum is non-degenerate the 
following inequality is true for any trial function i™^^{rLi, ^2) 

v^ix) ^ ^&(^^l, ^2)] < vo{x). (47) 

Choosing an appropriate trial function offers the possibility to determine a lower bound 
for vo{x) and therefore an upper bound for the critical length Ic. 

The boundary value problem described above is similar to those arising from 
electrodynamics or hydrodynamics. Its singular part implies that fnin2(^i)^2) is 
discontinuous at the contact surface. This is similar to the behavior of the electric 
field at a dipolar layer. Without restricting generality, we can choose R„2 — R^i along 
the 2;-axis. Then the relevant coordinate is = 02 — 0i and contacts can occur at 
= 0^ or = vr^ (see discussion in section 2), i.e. f„j„2(^ni, ^n2) = fnin2(^ij ^2, 0i, 0) 
must be discontinuous at = and = vr. Next, Eq. (|43| ) makes it obvious that 
|T,„^„2 ■ fnin2l should be made as large as possible. This can be done by choosing 
~ ui X = (cos 6*1 COS 01, COS 6^1 sin 01, — sin 6*1). Finally, the calculation 
ofjFfromEq. (^) requires the calculation of the r.h.s. ofEq. (|i5|). Since Vj/^^^^l^i; ^2) 
involves 1/ sin^j for x— and y— component of the gradient, the integrals in Eq. ( ^5] ) will 
not exist unless /"^^^ ^ (sin^i)^i(sin6'2)'^^ with z/i > 1 and 1/2 > 1. These considerations 
suggest to use the following trial function: 



f™:,^(fii,fi2) = \hZ^,{ni,n2) = Asin^isin^2 



X 



2 
2 



^ COS 6*1 COS 01 \ 

COS 6*1 sin 01 
1^ — sin 6*1 J 

1), O<02-01<7r 
bi), TT < 02 - 01 < 27r 



(4^ 



with A being a variational parameter. For the (02 — 0i)-dependence we have made the 
simplest choice taking a linear variation with discontinuities at and vr. Note also that 
fnin2(^i'^2) does uot depend on either R„2 — R„j or a; = //r„, which will be not true 
for the true, exact solution f^°'„^(fii, ^2)- Substituting f™5^^(fii, ^2) into Eq. (pf), where 
the first term is evaluated by Eq. (|45|), one obtains: 



^[f™^^(^]l,^]2)] = F{\) = \'h - 2A/i(x) 



with 



3(47r 



(49) 



(50) 
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1 1 

'3(47r)2 



X 



Note that I2 does not depend on x. F{X) is easily maximized. As a result one obtains: 

vrix) = F{X^^.^^) = ^-^^ . (52) 

I2 can be calculated analytically: 

/2 = -^(27 + 27r2) ^ -1.154 (53) 

ol 

whereas /i(a;) must be calculated by numerical integration, due to the nontrivial contact 
hypersurface. The result i;™''(x) following from this approach is shown in Fig. |^. It is 
interesting that Vq^^{x) reproduces the crossover behavior found for Vq'^^Ix) close to 
X = 2 from v™{x) ~ for x < 2 to t;™''(x) > for x > 2. Comparison of i)o"™(a;) with 
v™{x) confirms that i;™''(x) is a lower bound for t;™™(x), as it should be. Introducing 
vo{x) from both approaches into Eq. (^) leads to v{l) represented in Fig. |^ for an 
sc-lattice. From this figure we get: 

^num^3 45 ^ /™'^^5.68 : sc lattice. 
v{l) can be calculated for any periodic lattice. For the other cubic lattices we obtained 



;num ^ f 2.78 ^ f 4.84 bcc lattice 

\ 2.20 ' " \ 3.96 fee lattice. 

Since t;™''(0 < t;'^"'"(/) = v{l), it follows that Z™'^ > = l^. The critical length 
decreases when going from sc to bcc to fee lattice. This decrease is related to (i) the 
increase of the coordination number z from 6 to 8 to 12 and (ii) a decrease of the nearest 
neighbor distances ri from 1 to \/3/2 to \/2/2 (in units of the lattice constant a). Since 
an increase of z and a decrease of ri results in an increase of the steric hindrance, Ic must 
decrease. As the increase of steric hindrance is equivalent to an increase of collisions, i. e. 



of contacts, our result has some similarity to that found in ref. |^9|. These authors have 
used a probabilistic approach in order to derive a criterion for the mechanical stability 
of an off-lattice system of infinitely thin hard rods in its randomly closed-packed glassy 
state. It has been found that a mechanically stable amorphous phase occurs if the 
average contact point per rod becomes 5. However, such an arrangement of hard rods 
with d = is dynamically unstable, as discussed in the Introduction. In the simulation 
performed in ref. it was found /™ = 2.7 for an fee lattice. On the other hand ref. 



found = 4.5 for an sc lattice. Note that in ref. [301 the rods were fixed with 



one of their endpoints, and not with their centers. Our result for Z"*^™ and Z™'' agree 
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satisfactorily to the simulational result of ref. p9|. The difference between our result 



and that of ref. ||3^ may be due to the different way of fixing rods on the lattice that 
was mentioned above. 

The authors of ref. have also studied the /—dependence of the rotational 

diffusion constant D{1) defined by: 



D{1) 



lim In 



N 



N 



^(u„(t)-u„(0)) 



n=l 



(54) 



The variation of D over about two orders in magnitude follows a power law D{1) ~ 
{Ic-iy" with 7*^^ = 4.2. Since v{l) is analytical for / > 1 it follows that l-v{l) ~ Ic-l 
and therefore D{1) ~ (Z^ — 0- Hence our analytical theory yields 7 = 1 which differs 
significantly from 7*^^. This deviation may have two reasons. First, our result 7 = 1 
is mean-field-like (see, e.g., ref. [^). Note that in ref. going beyond mean field 
approximation it has been found 7 = 7/6 which is still close to one. Second, 
^MD ^ 4 2 is unusually high. Simulations and experiments of supercooled hquids usually 
yield 7 



2 consistent with most MCT-analyses ||T^, |T5|, |T6|, |T8|, jl 



4.2. Nonergodicity parameter and dynamics close to Ic 

The analysis in subsection 4.1 has proven that our theory predicts a dynamical glass 
transition at a critical length /c- There are two important questions remaining. First, 
does the nonergodicity parameter change at Ic in a continuous (type-A transition fl^) 
or discontinuous (type-B transition |]14|) way? Second, what is the time- or frequency- 
dependence of the 1-rod density p^^-* close to /c? The present subsection will give answers 
to both questions. 

Let us expand p^\^; z) and j^j^''(^^; z) into spherical harmonics: 



A 



where A = {j,m), j 
approximation (^) to finite z: 



0,1,2,--- 



-j < m < j. 



(56) 

Furthermore, let us extend 



(D(l]i, z))''^ ^ D{z)5''^ 5{ni\n2). (57) 
It is easy to prove then that the Laplace transform of Eqs. (|15D and (|2l|) leads to 

Uz) = f^l\z)/p^l\t = 0) = l/[z + J(J + l)D{z)], (58) 

where we skipped the index n. 

The nonergodicity parameter f\ is defined as follows: 



/a 



lini z(t)x{z) 



(59) 
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Taking into account approximation Eq. (|57| ) we get from Eq. (^) after operating 
with (1/4 vr) / dO^i f dVt2 on both sides the self-consistency equation for D{z): 

^ + J-v{C:l), C = z/D{z), (60) 



D{z) Do D{z) 

with 



{;(C;0 = ^E(Ton-[C-/:^'^1-'Ton). (61) 



Inverse Laplace transform of Eq. (^) yields 

1 
3 



t^(^;0 = ^E(Ton-e^''''*Ton). (62) 



n 



Note that v{t; I) is the time-dependent analogue of v{l) defined in Eq. (^). 

Eq. ([60| ) has a similar mathematical structure as the corresponding equation for a 
Lorentz gas obtained from a mode coupling approximation |^0|, Assuming that F{t) 
does not decay slower than t~^/'^ for t oo [52] it follows from the Tauberian theorems 
H that 

v{U) = v{Q\l)+v\Q-J)C + 0{e/^) (63) 



where 



i;(0;/) = v{l) 



^^'(0; = -\ E(Ton ■ (/:(')^)-'Ton) := -a{l) < 0. (64) 

-J n 

Substituting Eq. (^^^ into (§D|) and neglecting 0(C^/^) we get a quadratic 
equation for z/D{z), the physical solution of which is 

^/D{z) = \ W) + ^e\l)+z.U] (65) 

with 

(n - ^(^) ~ ^ _ [ < 0' ^ < (ergodic phase) 

~ a(/c) ~ 1 > 0, / > /e (glass phase) ^ ' 

and to = 4/(a(/c)-Do) being a microscopic time scale. Note that Eq. ( |65D is identical to 
the corresponding equation obtained for the so-called Fi-model []T^. Accordingly, all 
results for the Fi-model hold here as well. This implies (i) that the glass transition 
is of type A, i.e. the nonergodicity parameters (obtained from Eqs. (|58|) and (p9|)) 



, l=<l 



c 



fx{l) = { m ^ v'iQ , (67) 

vary continuously at Ic, (ii) at the critical point / = Ic and for z ^ 0, t oo 

M^) = hztoY^', Mt) = -^(t/to)-'/' (68) 

2 V47r 



and (iii) for |/ — /d <^ 1 and for z —>■ 0, t oo 

Mz, e) = (|e|M) 0a(^M), 0a (t, e) = MMt/Q (69) 
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with uj^ = = e^/to- In other words, the {z, e)- and (t, e)-dependence follows from the 
master function (j)x and (f)x (Eq. (|68D), respectively, by use of the scaled frequency z/lj^ 
and time t/t^. The time scale diverges as when approaching the glass transition 
independent on the sign of e. 



5. SUMMARY AND CONCLUSIONS 



Our main motivation has been the microscopic description of the glass transition for 
systems with trivial statics, i.e. systems which do not undergo an equilibrium phase 
transition to an ordered phase and do not have static correlations. For such systems, 
present microscopic theories like MCT and replica theory for structural glasses predict 
neither a dynamical nor a static glass transition. As a model we have chosen N infinitely 
thin hard rods with length L fixed with their centers on a periodic lattice with lattice 
constant a. The only relevant physical parameter is the dimensionless length / = L/a. 
Simulations for an fee ||2^ and an sc lattice |^ strongly suggested that a dynamical 
glass transition occurs at a critical length = 2.7 ^ and 11'"^ = 4.5 Note that 
the rods in ref. were fixed at the lattice sites with one of their end points. 

To describe the dynamics of our model we have used the generalized A^-rod 
Smoluchowski equation from which we derived a hierarchy of coupled equations for 
Vnr--nfe(^i' ' ' ' ' ^k), the fluctuations of the reduced fc-rod densities. Truncating at the 
second level and approximating the influence of surrounding rods on a pair of rods by 
the introduction of an effective diffusion tensor D we finally obtained a self-consistency 
equation for D (Eq. (pS])). This equation describes a feedback mechanism (as MCT also 
does) which leads to a slowing down of the dynamics and ultimately to an orientational 
glass transition, similar to that occuring in plastic crystals g |, 0, |, |ll], P, [H. The 



model studied in the present paper may be applied to real plastic crystals as has been 
shown in ref. p. 

One of the essential features of present MCT for the glass transition is the cage 
effect |T^. For low temperature or high density an extended particle, e.g. a hard sphere, 
is captured in a cage. The occurrence of this cage is accompanied by the growth of static 
correlations. Particularly, the main peak of the static structure factor grows and this 
leads to an increase of the static vertices entering the MCT-equations. If the vertices 
reach a critical strength the system will undergo a dynamical glass transition, provided 
ordering is prevented. In this sense the cage effect entering into MCT through the 
static correlations is of static nature. The model studied here does not have any static 
correlations, yet there is a cage effect. If the rods are long enough the tagged rod's 
motion can be strongly restricted to an "orientational" cage (cf. Figure 2). This cage, 
however, is of pure dynamical nature and leads to a glass transition as well. 

Comparing the results for the fee lattice obtained from our theoretical framework 
with the additional approximation Eq. (^Tj) with those from simulations 3D| one can 
say that both values = 2.20 and Z™'^ = 3.96 are in a reasonable range of = 2.7 
29| . Taking into account that our approach involves uncontrollable approximations 
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(Eqs. (|T9|) and (|20|)) this is a satisfactory result. The tendency of Ic to decrease 
when going from the sc lattice to the bcc to the fee lattice is consistent with our 
expectation that Ic decreases with increase of the coordination number and decrease 
of the nearest neighbor distances. Of course, second-, third-, etc. nearest neighbors 
play also a role, but with minor influence. Less satisfactory is the exponent, 7 of the 
power law, D{1) ~ {1^ — ly, which is 7 = 1, in contrast to 7*^'™ = 4.2 Even if 

the simulated result turns out too high, the type- A transition result, 7 = 1, cannot be 
the correct value. This requires an improvement of our microscopic approach. It is not 
obvious whether the approach presented in ref. which is based on rather qualitative 
arguments can serve as a guide. 

It seems that the "dynamical cage effect" leads to a discontinuous variation of the 
nonergodicity parameter at Ic and is accompanied close to Ic by a two-step relaxation 
process, as can be seen from the orientational correlator studied in ref. [|^. Such a 
two-step relaxation was first predicted by MOT [||, |1^, |15|, |T^ and within this theory 
it is related to the cage effect. The present theoretical approach predicts a continuous 
transition for the glass order parameters f\. Therefore it does not yield a two-step 
relaxation. This fact demands for an improvement of our theory, based on the hierarchy 
of equations for the k-rod densities. But on the other hand it is, of course, a big challenge 
to extend the present MOT such that it includes the "dynamical cage effect" described 
above and that a discontinuous glass transition occurs. If this turns out to be possible, it 
seems that the vertices must be generalized. Besides static correlations, i.e. correlators 
at time t = 0, they also must contain dynamical ones. These dynamical correlations 

may enter through a time- dependent force-force correlator / dt{F{t) ■ F(0)), i.e. as a 



Laplace transformed correlator at frequency z = (as found in the present approach). 
For systems with vanishing static correlations, like for our model, the vertices would 
be of purely dynamical nature. Now, increasing d to finite values will generate static 
correlations. For small thickness the "dynamical cage" will be still dominant, but at a 
crossover value d^o. the static cage effect would become comparable with the dynamical 
one and for d > d^.o. it would be the dominant one. 

Let us come back to the choice of Brownian vs. Newtonian dynamics. We have 
found that the glass transition is driven by the increase of the Laplace transform at z = 
of the time dependent torque-torque correlator of an isolated pair of rods. This correlator 
is be different for Brownian and Newtonian dynamics. Therefore the critical length Ic 
will be not the same for both dynamics. At present, it is not clear whether this is only a 
small effect. This is different from what has been found for liquids with nontrivial static 
correlations []33| , [3^ . However, the linear dependence of the nonergodicity parameters 
and of the diffusion constant on e, the quadratic dependence of the frequency scale 
cUe on e and the t"^/^ dependence of the correlator at the critical length Ic should be 
independent of the type of djTiamics. 

It would be interesting to re-investigate the lattice model with infinitely thin hard 
rods by computer simulations. With present computers it should be possible to cover 
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a larger time range and to study the dynamics close to Ic in greater detail. This would 
allow to check whether the two scaling laws and other predictions of MCT are consistent 
with simulations. If the outcome of such a test proves consistency it would encourage 
theoretical efforts to extend the microscopic theory within the framework of MCT. 

Besides rods on a lattice one could also investigate liquid systems of "particles" 
built by crossing infinitely thin hard rods. If the number of "legs" (rods) of a "particle" 
becomes large they may exhibit quite similar dynamical behavior to hard spheres, 
although there are no static correlations. The time scale on which these "particles" 
realize that they are not hard spheres could be much larger than a typical time scale 
for structural relaxation. 

To conclude, we have discussed a purely dynamical mechanism which also drives 
a glass transition and which up to now has received almost no attention. Its further 
investigation by computer simulations and analytical work seems to us a challenging 
task for the future. 
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Figure 1. Unit cell of the sc-lattice with lattice constant a with hard rods of length 
L. 



Figure 2. Illustration of the blocking of rod 1 by rods 2-5 within shaded region for 
/ > 2.. 



Figure 3. Illustration of possible collisions of rod i and j for fixed rod i. For better 
visualization a small but finite thickness has been used, (a) Single collision for Z < 2 
at = 4>f . The rod j' (dashed) turned around z-axis by tt does not lead to a second 

collision if / < 2. (b) Single collision for Z < 2 at = (f>f + tt. Again a rotation by tt 
(not shown) does not yield a second collision (c) Two collisions for Z > 2; 0j and 0j are 
chosen such that two collisions occur at = (j)f and = (j)f + 77. 



Figure 4. Illustration of geometrical quantities defined in the text for two rods i and 
j- 



Figure 5. Time-dependence of the torque-torque correlator of an isolated pair of rods 
for different Z/r„. Dotted line: Z/r„ = 1.8; dash-dotted line: Z/r„ = 2; dashed line: 
Z/r„ = 4; solid line: Z/r„ = 6 



Figure 6. vo{x) from the numerically exact (squares) and the variational calculation 
(solid line). 



Figure 7. v{l) for a sc-lattice from the numerically exact (squares) and the variational 
calculation (solid line). Zc denotes the critical length for which u(Zc) = 1. 



